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^ ', Abstract 

The purpose of this paper is to find optimal estimates for the Green function of a 
half-space of the relativistic a-stable process with parameter m on M. d space. This process 
has an infinitesimal generator of the form ml — [m 2 / a I — A) a / 2 , where < a < 2, m > 0, 
^ ' and reduces to the isotropic a-stable process for m = . Its potential theory for open 

CN ■ bounded sets has been well developed throughout the recent years however almost nothing 

was known about the behaviour of the process on unbounded sets. The present paper is 
intended to fill this gap and we provide two-sided sharp estimates for the Green function 
Ph ■ for a half-space. As a byproduct we obtain some improvements of the estimates known 

^ for bounded sets. 

Our approach combines the recent results obtained in [5], where an explicit integral 
formula for the m-resolvent of a half-space was found, with estimates of the transition 
densities for the killed process on exiting a half-space. The main result states that the 
Green function is comparable with the Green function for the Brownian motion if the 
' points are away from the boundary of a half-space and their distance is greater than one. 

On the other hand for the remaining points the Green function is somehow related the 
Green function for the isotropic a-stable process. For example, for d > 3, it is comparable 
with the Green function for the isotropic a-stable process, provided that the points are 
\Q ■ close enough. 
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1 Introduction 

In the paper we deal with some aspects of the potential theory of the a-stable relativistic process. 
That is a Levy process on M. d with a generator of the form 

H™ = ml - (m 2/a I - A) o/2 , < a< 2, m > 0. 

For m = the operator above reduces to the generator of the a-stable rotation invariant 
(isotropic) Levy process which potential theory was intensively studied in the literature. 
For a = 1 the operator 

H™ = ml- (m 2 I - A) 1 / 2 
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plays a very important role in relativistic quantum mechanics since it corresponds to the kinetic 
energy of a relativistic particle with mass m. Generators of this kind were investigated for 
example by E. Lieb [19] in connection with the problem of stability of relativistic matter. An 
interested reader will find references on this subject e.g. in a recent paper [T8] . 

Another reason that the operator H™ is an interesting object of study is its role in the 
theory of the so-called interpolation spaces of Bessel potentials and its application in harmonic 
analysis and partial differential equations (see, e.g. [23] and p3]). This theory is based on 
Bessel potentials defined as J a — (I — A)~ Q / 2 . As Stein pointed out in his monograph [23], the 
Bessel potentials exhibit the same local behaviour (as \x\ —>■ 0) as the Riesz potentials but the 
global one (as \x\ — > oo) of J a is much more regular. In terms of the relativistic process the 
potential J a is so-called 1-resolvent kernel of the semigroup generated by H^. 

In the paper we consider the process killed on exiting the half-space H = {x G R d : x d > 0} 
and examine the behaviour of its Green function Gu{x, y)- Contrary to the stable case a closed 
formula for that Green function is not know and seems to be a very challenging target. Recently 
in [5J an integral formula in terms of the Macdonald functions was found for G^(x,y) - the 
m-resolvent kernel for H. As proved in [3], for d > 3, the behaviour of the Green function is 
equivalent to the behaviour the m-resolvent if \x — y\ — > 0. Our main result establishes optimal 
bounds for the Green function of HI. To our best knowledge it is the first result of that type 
when optimal estimates for unbounded set (different than the whole M. d ) are derived. 

At this point let us mention that the potential theory for bounded sets has been well 
developed during recent years (see [S], [22], [H], [ED- Under various assumptions of the 
regularity of a bounded open set D it was shown that the Green function of D was comparable 
with its stable counterpart. This comparison allowed to prove the relativistic potential theory 
shares most of the properties of the stable one if bounded sets are considered. Comparing the 
potential kernel for the stable process with the potential kernel for the relativistic process (see 
[20J) we can conclude that such a comparison of Green functions is not generally possible for 
unbounded sets. Since the relativistic potential kernel (for d > 3) is asymptotically equivalent 
(if \x — y\ is large) to that of the Brownian motion it may suggest that the Green function of H, 
at least for some part of the range of x, y, is comparable with the Green function of H for the 
Brownian motion. Our main result confirms that suggestion and we prove the comparability 
for points x, y being away from the boundary and with \x — y\ > 1. For other points our bound 
is also optimal. 

We also thoroughly examine the one-dimensional case and provide optimal estimates for the 
Green functions for bounded intervals taking into account their length. While for intervals of 
moderate length (say smaller than 1) we can use the well known results about comparability 
of stable and relativistic Green functions, for large intervals we relay on the estimates for half- 
lines obtained in this paper. Again we show that the Green functions for large intervals are 
comparable to the Brownian Green functions for most of the range. 

The organization of the paper is as follows. In Section 2 we collect all definitions and prelim- 
inary results needed for the rest of the paper. The next section is basic for the paper. Here we 
prove the estimates for the Green function of (0, oo). Then in Section 4 we apply them to prove 
the optimal bounds for the tail function of the exit time from (0, oo) and some other properties 
of the exit times. These estimates will have a crucial role in examining multidimensional case 
which was accomplished in Section 5. We conclude the paper with exploring in the last section 
the one-dimensional case with regard to optimal estimates for bounded intervals. 
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2 Preliminaries 



Throughout the paper by c,C,C\ ... we denote nonnegative constants which may depend on 
other constant parameters only. The value of c or C, C\ ... may change from line to line in a 
chain of estimates. 

The notion p(u) ~ q(u), u G A means that the ratio p(u)/q(u), u G A is bounded from 
below and above by positive constants which may depend on other constant parameters only 
but does not depend on the set A. 

We present in this section some basic material regarding the a-stable relativistic process. 
For more detailed information, see [22] and [7]. For questions regarding Markov and strong 
Markov property, semigroup properties, Schrodinger operators and basic potential theory, the 
reader is referred to [9] and [3]. 

We first introduce an appropriate class of subordinating processes. Let 9 a (t,u), u, t > 0, 
denote the density function of the strictly a/2-stable positive standard subordinator, < a < 2, 
with the Laplace transform e~ txa/2 . 

Now for m > we define another subordinating process T a (t, m) modifying the correspond- 
ing probability density function in the following way: 

9 a {t, u, m) = e mt 6 a {t, u) e ~ m2/au , u>0. 

We derive the Laplace transform of T a (t,m) as follows: 

E e -XT a (t,m) = e mt e -t(A+m*/<*r/^ A > _ m 2/a (2 ^ 

Let B t be the symmetric Brownian motion in IR d with the characteristic function of the form 

E o e it-B t = e -t\e (2.2) 

Assume that the processes T a (t,m) and B t are stochastically independent. Then the process 
Xt' m = BT a (t,m) is called the a-stable relativistic process (with parameter m). In the sequel we 
use the generic notation X™ instead of X" ,m . If m = 1 we write T a (t) instead of T a (t, m) and 
X t instead of X\. From (12.11) and (I2.2p it is clear that the characteristic function of X™ is of 
the form 

Obviously in the case m = the corresponding process is the standard (rotationally invariant 
or isotropic) a-stable process. X 4 m is a Levy process (i.e. homogeneous, with independent 
increments). We always assume that sample paths of the process X 4 m are right-continuous and 
have left-hand limits ("cadlag"). Then X t m is Markov and has the strong Markov property 
under the so-called standard filtration. 

From the form of the Fourier transform we have the following scaling property: 

p?(x)=m d / a pU™ 1/a x). (2.3) 

In terms of one- dimensional distributions of the relativistic process (starting from the point 0) 
we obtain 

X? ~ m-^Xmt , 

where X t denotes the relativistic a-stable process with parameter m = 1 and "~" denotes 
equality of distributions. Because of this scaling property, we usually restrict our attention to 
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the case when m — 1, if not specified otherwise. When m = 1 we omit the superscript "1", i.e. 
we write pt(x) instead of p\ (x), etc. 

Various potential-theoretic objects in the theory of the process X t are expressed in terms 
of modified Bessel functions K v of the second kind, called also Macdonald functions. For 
convenience of the reader we collect here basic information about these functions. 

K u , v G M, the modified Bessel function of the second kind with index z/, is given by the 
following formula: 

r°° 2 

K„(r) = 2- 1 - u r u / e^e'^u-^du , r > 0. 
Jo 

For properties of K v we refer the reader to [TT]. In the sequel we will use the asymptotic 
behaviour of K v : 



K v (r) = ^C-) . ; - 0+. //-(). (2.4) 



2 V2, 

K (r) = -logr, r^0 + , (2.5) 



K v {r) = ^Le~\ r^oo, (2.6) 
y2r 

where g(r) = f(r) denotes that the ratio of g and / tends to 1. For v < we have K u (r) = 
K_ u (r), which determines the asymptotic behaviour for negative indices. 

The a-stable relativistic density (with parameter m = 1) can now be computed in the 
following way: 

POO 

Pt(x)= / e t e a {t,u)e~ u g u {x)du, (2.7) 
Jo 

where g u (x) = j^~jd/2 e ~^~ is the Brownian semigroup, defined by (12. 2p . 

We also recall the form of the density function u(x) of the Levy measure of the relativistic 
a-stable process (see [2"2"]): 

poo 

= 2 r(l - i) Jo e ~ U 9u(x)u- l - a/2 du (2.8) 

Oil 2 d+q , , , 

= 7r^r(i-f) |xr ^^ (N) - (2 ' 9) 

In the case < a < 2 we have the following useful estimates (see [22] for the proof of the 
first lemma): 

Lemma 2.1. There exists a constant c = c(a,d) such that 

maxp t (x) < c{t- d/2 + r d,a ) . (2.10) 

xeR d 

Lemma 2.2. For any t > and x G M. d we have 

p t (x) < c(d,a) (g t (x/V2) + tu(x /y/2] 

and 

, . c(d) 

Pt{x) - W' 
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Proof. Notice that for u, t > 0, 

e a (t, u) = r 2/a e a {i, r 2/a u) and e a (i, u) < cu^- a/2 . 

Hence 

6 a (t, u) < ctu- 1 -^ 2 , t,u>0. (2.11) 
Using (12. 7p we obtain for t > 1, 

2 /-2t 

p t (x) < e-^(47r)-^e* / e a (t,u)e- u u- d/2 du 



+ ct g u {x)e 2« 1 a/2 <iu 

2 , /"OO 

< e - « (47r)-3e* / 6 a (t,u)e- u U - d/2 du 



+ ct g u {x)e au q/ 

= e 8' pt(0) + cbv(x/y/2), 
where we used ( 12. 8 p in the last line. Moreover by Lemma 12.11 we can estimate 

e-^pt(O) < cg t (x/V2), t > 1. 



This completes the proof of the first estimate for t > 1. 

Next, for t < 1, applying fpTT]) . (I2~TTT) and (EHD we arrive at 

/■oo 

p t (x) < c£ / ^ n (x)e""M" 1_Q/2 rfM = ctz/(x) < ctu(x/V2), 



which complete the proof the first inequality. 

The second bound is true for the transition density of any subordinated Brownian motion. 
Indeed let us observe that for any t > and x G M. d , 

, . / d \ d ^ 2 _d . , , 
gt{x) < r- e a ar d . 



2tt y 

Hence by subordination 

ef/2 



p t (x) = Eg Ta{t) (x) < (j^J e t\x\ 



□ 



The standard reference book on general potential theory is the monograph [3J. For conve- 
nience of the reader we collect here the basic information with emphasis on what is known (and 
needed further on) about the a-stable relativistic process. 

In general potential theory a very important role is played by A-resolvent (potential) kernels, 
A > , which are defined as 



oo 

At_, /„ „,\sJ+ ™ „. r- inxi 



U\(x,y)= I e p t (x-y)dt, x,y e 
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If the defining integral above is finite for A = 0, the corresponding kernel is called a potential 
kernel and will be denoted by U(x,y). For the relativistic process the potential kernel is well 
defined for d > 3 but contrary to the stable or Brownian case it is not expressible as an 
elementary function. Recall that for the isotropic a-stable process the potential kernel is equal 
to C\x — y\ a ~ d for d > a and for the Brownian motion it is C\x — y\ 2 ~ d for d > 3, where C's are 
appropriate constants. One can prove that the relativistic potential kernel could be written as 
a series involving the Macdonald functions of different orders but this formula does not seem 
very useful. Nevertheless the asymptotic behaviour of the potential kernel was established in 

U{x — y) « \x — y\ a ~ d , \x — y\ < 1; U(x - y) « \x - y\ 2 ~ d , \x - y\ > 1. (2.12) 

Note that they suggest that the process locally behaves like a stable one and globally like a 
Brownian motion. Despite the fact we do not know any simple form for the potential kernel, a 
formula for the 1-potential kernel is known (e.g. see [5]): 

TT ( \ nt A\ K (d-a)/2{\x\) 

U 1 (x) = C(a,d) |a;[(d _ a)/2 , 

where C{a,d) = r(a/2) ^ /2 • 

The first exit time of an (open) set D C M. d by the process X t is defined by the formula 

r D = M{t > 0; X t i D} . 

The basic object in potential theory of X t is the X-harmonic measure of the set D. It is 
defined by the formula: 

P£(x,A) = E x [t d < oo;e- x ^l A (X TD )]. 

The density kernel of the measure P^(x,A) (if it exists) is called the X-Poisson kernel of the 
set D. If A = the corresponding kernel will be denoted by Pd{x,z) called Poisson kernel of 
the set D. 

Another fundamental object of potential theory is the killed process Xf when exiting the set 
D. It is defined in terms of sample paths up to time tjj. More precisely, we have the following 
"change of variables" formula: 

E*f(X t D )=E*[t<T D ;f(X t )], t>0. 

The density function of transition probability of the process Xf is denoted by . We have 

pf (X, y) = Pt (x -y)-E x [t> r D ;p t _ TD (X TD - y)\ , x,yeR d . (2.13) 

Obviously, we obtain 

Pt(x,y)<p t (x,y), x,yeR d . 

pf is a strongly contractive semigroup (under composition) and shares most of properties 
of the semigroup p t . In particular, it is strongly Feller and symmetric: p^(x,y) = p^(y,x). 

The A-potential of the process Xf is called the X- Green function and is denoted by G^. 
Thus, we have 



i>oo 

G x D (x,y)= / e- xt p?(x,y)dt 
Jo 
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If A = the corresponding kernel will be called Green function of the set D and denoted 
G D (x,y). 

Integrating (12.131) we obtain for A > 0, 

G x D (x, y) = U x (x, y) - E x e- Xr °U x (X TD ,y). 

Suppose that D\ C D 2 are two open sets. By the Strong Markov Property 

G x D2 (x,y)-G x Dl (x,y) 
= E* [e- XT "iU x (X TDi ,y) - e^U x (X TD2 ,y 

= E 2 
= E 2 



r Dl < r D2 ;e- x ^ (u x (X TDi , y) - e~ XTD ^ U X (X TD2 , y) 
r Dl < r D2 ;e- Xr ^ (u x (X TDi ,y) - E X ^e~ Xr ^U x (X TD2 ,y) 



= E x [r Dl < T D2 ;e- x ^G x D2 (X TDi ,y)] . (2.14) 

The main purpose of the present paper is to obtain sharp estimates for the Green function 
for D = M = {x G M. d : x^ > 0}. The investigation of Green functions of the relativistic 
process for unbounded sets seems not to be treated in the literature. For bounded sets there 
many results obtained in recent years showing that the Green functions for open bounded sets 
under some assumptions about regularity of their boundary are comparable to their stable 
counterparts in M. d , d > a ([22], [H], [TTj ) . That is, for x,y G D, 

CiD^Gg^fay) < G D (x,y) < C(D)G% able (x,y), (2.15) 

where G s ^ able is the corresponding Green function for the isotropic stable process and C(D) is a 
constant usually dependent on diam(D) = sup Xty£ £i\x — y\. Unfortunately in all known general 
bounds of the above type the dependence on the set D in the constant C(D) is not very clear 
and C(D) grows to 00 with diam(D). The constant also depends on some other characteristics 
of D as e.g. Lipschitz characteristic of D when D is a Lipschitz set. Therefore it is not possible 
to use well known exact formulas or estimates for the stable Green functions of regular sets 
as half-spaces, balls or cones to derive the corresponding optimal estimates for the relativistic 
process. Even for balls the constants grow to 00 and ( 12. 15ft does not yield any estimate for a 
half-space in the limiting procedure. 

Now suppose that D is a bounded set with a C 1 ' 1 boundary. It is well known that there is 
a p > such that for each point z G dD there are balls B z C D, B* C D c of radius p such 
that z G B z fl B*. Denote by p = po(D) the largest p having the above property. Finally let 
7 = diamD/po- However not explicitly stated, the following bound can be deduced from the 
results proved in [22], for d > a, x,y G D: 

C 1 ( 7 )C(diam( J D))- 1 Gg aMe (x, Z /) < G D (x,y) 

< C 2 ( 1 )C(di a m(D))G^ able (x,y), 

(2.16) 

where the constant C can be chosen in such a way that C(diam(D)) = 1 for diam(D) < 1 and 
C(diam(Z))) increases with diam(D). With some extra effort one can prove that the growth is 
polynomial. The constants 6*1(7), 6*2(7) can be chosen as continuous with respect to 7. Note 
that if D is a ball than we can take absolute constants (depending only on a and d) instead of 
61 (7), 62(7). 
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Hence for "smooth" sets with small or moderate diameter the estimate (I2.16P is very satis- 
factory. For example for balls of small or moderate diameter we obtain very precise estimates 
using well known results for the isotropic stable process. However, in the case of balls of large 
size, it would be very interesting to find optimal estimates of the relativistic Green function. 
Our main result provides optimal estimates for the Green function of the half-space H. Also we 
found optimal estimates for intervals in K. Despite the fact we do not examine Green functions 
for balls in higher dimensional spaces we provide very precise estimates of the expected first 
exit time from a ball. 

Now we define harmonic and regular harmonic functions. Let u be a Borel measurable 
function on K d . We say that u is harmonic function in an open set D C R d if 

u{x) = E x u(X TB ), x G B, 

for every bounded open set B with the closure B C D. We say that u is regular harmonic if 

u(x) = E x [t d < oo; u(X TD ))], x e D. 

As a result of (I2.16P we obtain the following version of the Boundary Harnack Principle (for 
details see [22] or [H] in the one-dimensional case). 

Theorem 2.3. [BHP] Let D be a bounded set with a C 1,1 boundary. Suppose that diamD < 4 
and po(D) > 1. Let z G 3D. If f is a non-negative regular harmonic function on D and 
fix) =0,i6 B(z, 1) n D c . Then 

f(x) « f(x )6 D (x) a / 2 , xeB(z,l/2), 

where 5rj{x) = dist(x, dD) and x G D such that 5d{xq) = 1. 

For the purpose of this paper we state the following specialized form of BHP which can be 
easily deduced from Theorem 12.31 

Lemma 2.4. Let H 3 1 = (0, . . . , 0, 1) and let F = B(0, y/2) n EL Suppose that f is a regular 
nonnegative harmonic on F such that f(x) = 0, x G HI C . Then for every x G B(0, 1) fl EI we 
have 

f{x) » fil)xf. 

Assume that R>2. Let D = 5(0, R), z = (0, . . . , 0, R) and x = (0, . . . , 0, R - 1). Let F = 
.B(zO) 2)(~\D. Suppose that f is regular nonnegative harmonic on F such that fix) = 0, x G D c . 
Then for every x G -B(^o, 1) H D we have 

fix)^fix )iR-\x\) a ' 2 . 

As mentioned above, the one-dimensional case for intervals was treated recently in [14] and 
since we will need it in the next section we present it in a convenient form of the estimate of the 
Poisson kernel. Actually in [T3j it was shown that the Green function of (0, R) is comparable 
with the Green function of the corresponding stable process (with uniform constant for R < 3). 
By standard arguments (see [22]) this implies the lemma below. 

Lemma 2.5. Assume that d = 1 and < R < 3. Let D = (0, R). Then 

P D ix, z) ps , r ,} X ^ R r ,., X }), r e~ z , x G D, z > R. 

v 1 (R(z - R)) a / 2 iz - x) 
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This implies that 

P X (X TD > R) « (x/R) a/2 , xeD 

and 

E X [X TD > R- X TD ] « (x/R) a/2 ((R - x) a/2 +x), xE D. 

We also have that 

E x r D ^ (x(R-x)) a/2 , xeD. 

Obtaining any exact formulas for the Green function or the Poisson kernel even for regular 
sets seems to be a very hard task but in the recent paper [S] the formulas for the 1-Poisson and 
1-Green function of H were described explicitly in terms of the Macdonald functions: 

Theorem 2.6. Let 

E x [e- m ,X m edu] = P^(x,u) 
be the 1-Poisson kernel for H. Then we have 

P&x,«) = 2 Sin(mr/ V)-^ ( ^X 12 Ki '^ X T^\ 
^ \—Ud J \x — u\ d l l 

•where Ud < < Xd- Let G^(x, y) be the 1-Green function for HI then for x, y e HI, 

G m (x,y)= (27r)rf/2r(a/2)2 J o ^^Kd/^-yKt + iY^dt. 

Moreover, 

[ G 1 m (x,y)dy = 1-E x e 
Jm 

r O/2) Jo 

This result will be very useful in our analysis since, as shown in [5] the behaviour of the 
Green function Gm(x,y) could be described in terms of the 1-Green function G^(x, y) when x 
and y are close enough. 

One of our main tools in establishing the upper bounds of the Green function will be 
estimates for the tail function P x {t^ > t). We start with the following lemma taken from the 
Master Thesis of the first author [13J. 

Lemma 2.7. There is a constant C such that 

P x (m>t)<c Xd + l ^ t 2 +1) , t>l,x d >0. (2.18) 

Proof. Let Y t = x[ d \ where X t = {X[ , . . . ,x[ ). By the symmetry of the random variable 
Yt we obtain 

P x (m>t) = P x (MY s >0) 

s<t 



-m 



P°(inf(-y s + x d ) > 0) = P°(sup Y s < x d ). 

s <t s<t 
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Using a version of the Levy inequality ([2J, Ch.7, 37.9) we have for any e, y > that 

n 

2P°{Y t >y + 2e) - 2V P°{Ytk - Y H k-p > e) < P°(su P r« > y). 

n » k<n n 

k=l — 

Note that Y^k=i P°(Ytk ~ Y t (k-i) > e) = nP°(Y± > s) —> t f°° is(x)dx, hence, by symmetry 



•■oo 

)£ 

•"71 

again 



POO 

P°(supF s > y) > 2P°(Y t >y + 2e)-2t / u(x)dx 

s<t Je 

/oo 
v(x)dx. 

This implies that 

POO 

P x (r m >t) = P°(sup Y s < x d ) < P°(\Y t \ <x d + 2e) + 2t v{x)dx . 

s<t Je 

For e > 1 we obtain from (I2.9P and (12.61) 

/oo 
u(x)dx < Ce- £ e- a/2 -\ 

Lemma |2~T1 implies that the density of Y{t) is bounded by Ci -1 / 2 , t > 1, hence taking e = 
| ln(i + 1) we obtain 

P x (re>t)<C(x d + ln(t + l))r 1 / 2 . 

□ 

In order to improve the above estimate for x close to the boundary we use Lemma 12.51 
Lemma 2.8. For < x d < 2 we have 

P x {m >t)< Cx a J 2 ln(t + I)/* 1 / 2 , t > 2, (2.19) 

where C is a constant. 

Proof. It is enough to prove the claim for d = 1. Let D = (0, 2) and assume that < x < 2. 
By the Strong Markov Property and then by Lemma [2.71 we obtain for t > 1: 

P x (r ( o,oo) > 2t) < P x {t d > t, r (0 ,oo) > 2t) 

+ S x [td < r (0)00 ) ; P Xt ° (r (0lOO ) > *)] 

< P x -(r D >t) 

+ ^[r fl < r (0 ,oo);X TD +ln(t + l)]/t 1/2 

< ^ + CP-[X TD >2;X TD ]/t 1 / 2 

+ CTn(t + l) P X '(X TD > 2)/t 1 /2 

< Cx a/2 \n(t+ l)/t 1/2 . 

The last inequality follows from Lemma 12.51 The proof is complete. □ 



10 



The estimates from Lemmas 12.71 12.81 will be very useful for the estimates of the Green 
function of the half-line, however they are not optimal. In the sequel we will be able to improve 
them to be sharp enough and optimal (see Proposition 14. 5ft . This will have a great importance 
in estimating the Green function for a half-space in the d- dimensional case. 

Lemma 2.9. There is a constant C such that for any open set D: 

Eg° a{t) (x,y) < p?(x,y) < C(t~ d l 2 + r d / a )P%r D > t/3)P^r D > t/3), 

where gf(x,y) is the transition probability for the Brownian motion killed on exiting D. 

Proof. We start with the upper bound. Since pf{x, y) is a density of a semigroup and pj?(x, y) < 
max 06K j Pt{z) then we have 

P3t( x iV) = P?(x,z)p^(z,w)p^(w,y)dzdw 
Jd Jd 

D/„ \ r / D 



< maxpt(z) / p t (x,z)dz / p t (w,y)dw 

2GM d J D J D 

= maxp t (z)P x (T D > t)P y (r D > t), 



which proves the upper bound since max zeR d p t (z) <C(t d l 2 + t by Lemma [2.11 

To get the lower bound we use the subordination of the process to the Brownian motion: 
X t = B Tai t). Then 

p°{x,y) = P x (B Ta{t) edy,B Ta{s] eD,0<s<t) 
> P x (B Ta(t) edy,B s eD,0<s< T a (t)). 

Using the independence of T a and the Brownian motion B we obtain 

P x {B Ta(t) e dy,B s eD,0<s< T a (t)\T a (-))=g° a[t) (x,y), 

Integrating we obtain the lower bound. □ 

The following lemma provides a very useful lower bound. Its proof closely follows the 
approach used in [20J, where the bounds on the potential kernels (Green functions for the 
whole M d ) were established for some special subordinated Brownian motions (in particular for 
our process) for d > 3. 

Lemma 2.10. For any open set D 6 M. d we have 

2 

G D (x,y) > -G 9 £ uss (x,y), 
a 

where G g D auss (x,y) is the Green function of D for the Brownian motion. 

Proof. Let Q(x,y) = Eg® , t Jx,y)dt. From the previous lemma it is enough to prove that 
Q(x,y)>lG 9 D auss (x,y). We have 

POD 

Q{x,y) = / Eg^ a{t) (x,y)dt 
Jo 

i POO 

e l I g®(x,y)e~' u 6 a (t,u)dudt 
Jo 

i POO 

g®(x,y)e~ u / e l 9 a (t,u)dtdu 
Jo 

i 

n 



where G(u) = e~ u J °° e t 9 Q {t, u)dt is the potential kernel of the subordinator T a (t). It was proved 
in [20] that G(u) is a completely monotone (hence decreasing) function and ini u>0 G(u) = 
lim^oo G{u) = C a . We find the constant C a by taking into account the asymptotics of the 
Laplace transform of G(u) at the origin: 



OO POO 



- Xu G(u)du = I e* J e- u{1+x) 9 a (t,u)dudt 
,t e -(i+\r/H dt 



(1 + A)«/ 2 



Applying the monotone density theorem we obtain that C a = 2/ a. Thus, since g®(x,y) > 0, 
we finally obtain 

POO 

Q{x,y) = / 9u( x ^y) G ( u ) du 



2 f°° 2 
> -/ 9u( x ,y) du = -G 9 D USS (x,y). 
a In a 



□ 



At this point let us recall that the exact formulas for the Brownian Green functions are well 
known for several regular sets as balls or half-spaces (see e.g. p]). Since some of them will be 
useful in the sequel we will list them for the future reference. Recall that the Brownian motion 
we refer to in this paper has its clock running twice faster then the usual Brownian motion. 
For the half-space H, for d > 3, we have that 



GZ ss (x } y) = C(d) 



1 1 



\x — y\ d ~ 2 \x — y*\ d ~ 2 



\x — y^' \x — y\ d ~ 2 



m H^^l i5=sh x,yem, (2.20) 



where y* = (y l , . . . , y d _ lt -y d ) G W. 
For the half-space H, for d = 2, 



Gr^) = 7flMl + 4^^), x,yeW. (2.21) 



27r \ |x — y 
In the one dimensional case 

G^'(x,y)=xAy, x,y>0. (2.22) 
For the finite interval (0, R) we have 

G%%(*,V) = X{R ~ V) R V{R ~ X) , x,ye{0,R). (2.23) 

Lemma 2.11. Let D be an open subset o/H. For fixed y £M the function Gu{-,y) is regular 
harmonic on D provided y ^ D. The same conclusion holds ifM is replaced by an open bounded 
set. 
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Proof. The proof is standard and is included for completeness. First observe that G-g.(z, w) < oo 
for z 7^ w, which follows from Lemmas 12.21 12.71 and 12.91 Next, applying (12.141) with D2 = H 
and D\ = D we have 

G x (x,y)-G x D (x,y) = E x [r D < r m ; e^G^^y)} . (2.24) 

If y D then t d = and X TD = y, P y - a.s. so Gp(x, y) = G^y^x) = 0. Moreover 
E x G^(X m ,y) = 0, which follows from the fact that P x (X m G H c \ (H c ) r ) = 0, where (H c ) r is 
a set of regular points of H c and for every z G (H c ) r , y G W 1 we have Gj^z, y) = (see [3]). 
This implies that (12.241) can be rewritten as 

G x l {x,y) = E x e~ x ^G x l {X TD ,y). 

Passing with A — > and observing that G^ /* G H we obtain the conclusion by the monotone 
convergence theorem. 

The same arguments can be applied for any bounded set F, since there is a half-space 
containing F, which guarantees that the Green function Gp{x,y) < 00 for x 7^ y. 

□ 

The estimates below following from Theorem 12.61 were proved in [BJ. They turn out to be 
useful in the next sections. 

Theorem 2.12. Assume that d = 1 and a > 1. When \x — y\>lAxAy>0we obtain 

Glo t0o) {x, y) « |, n l /9 (l AxA y)"/ 2 , 
u>/w/e /or \x — y\ < 1 /\ x /\ y we obtain 

G (0,oo)(^2/) ~ log 

G\ 0tOo) (x,y) « (lAxAy) 
In t/ie remaining case, a < d, we have 



x - yl 1 -*/ 2 

1 A x Ay 
\x - y\ 

a-l 



, if a = l, 
if a > 1. 



lf(«i-a)/2(|z-y|) 



F - 1/1 



i-a)/2 



1 AXg Ayg 

\x — y\ A 1 



a/2 



A 1 



Finally we state some basic scaling properties both for the Poisson kernel and the Green 
function. The proof employs the scaling property (12.31) and consists of elementary but tedious 
calculation hence is omitted. 

Lemma 2.13 (Scaling Property). Let D be an open subset o/M d and Pn,m, Gr,,m be the Poisson 
kernel, or the Green function, respectively, for D for the process with parameter m. Then 

P D , m ( x i u ) = m d/a P m i /aD (m 1/a x, m 1/a u), xeD,ueD c , 
G D , m {x, y) = m^l a G m , iaD [m x l a x, m^y), x G D,y G D . 
Thus, if D is a cone with vertex at we obtain: 

P D , m {x,u) = m d/a P D (m 1/a x,m 1/a u), x G D,u G D c , 
G D , m (x, y) = m^' a G D {w}/ a x, m^y), x G D,y G D . 
Due to these scaling properties it is enough to investigate the case m = 1. 
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3 Green function of half-line 

In this section d — 1 and the half-space EI is a half-line, that is H = (0, oo). 
Lemma 3.1. Assume that \x — y\ < 3. Then there is C = C(a) such that 

(lAxAy) a/2 <CGl 0>oo) (x,y). 

Proof. We use Theorem 12.121 First let a > 1 and \x — y\>lAxAy>0, then 

Glo l0o) (x,y) > ce-^x-y^ilAxAyT' 2 
> ce^T^-^lAxAyT 12 . 

Suppose that \x — y\ < 1 A x A y. For a — 1, 

1 A x Ay' 



Glo >oo) (x,y) > clog 



\x - y\ 



> clog 2 > c(l A x A y) 



a/2 



For a > 1. 



Gjo.co)^ y) > c(l A x A y) a ~ l > c(l A x A y) a ' 2 . 
Next, observe that Kn a \/2(r) /r^-~ a ^ 2 is decreasing. Therefore for a < 1 we obtain 



G\ Qi00) {x,y) > c 



K {1 _ a)/2 (\x - y\) 



| x _ y |(l-a)/2 



1 A A \ a /2 

1 Ax Ay \ 



| x — y\ A 1 



A 1 



> c- 



•^(l-a)/2(3) , , a / 2 

(1 A x A y) ' . 



3(l-a)/2 



□ 

Theorem 3.2. For x, y > 0, 

G {0>oo) (x, y) w G; 0iOo) (x, y) + (x A y) V (x A y) a/2 . 

Proof. Throughout the whole proof we assume that < x < y. The proof will rely on the 
estimates of P x (t(q }00 ) > t) derived in the previous section and the application of Lemma [231 
We proceed to estimate the Green function from above. First we split the integration 



Pi 



(0,00) 



x, y)dt 



pi°' 00 \x,y)dt + 
V(x,y) + R(x,y). 



Pt 



(0,00) 



x, y)dt 



We start with the estimation of the second integral. Due to Lemma 12.91 



R(x,y) 



P3t°°\ x iy) dt 



r°° dt 
< c / P x (r { o !Qo) > t)P ?/ (^(0,oo) > t) — . 

J 2 

First consider the case y < \A2- Then using (12.191) we have 



R{x,y)<C{xy) a ' 2 



(hit) 2 dt 



t t 1 / 



a/2 
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If y > v2, using (12.181) we estimate 

/■OO 7 . 

R&y) < c (P y (r i0 , oo) >t)) 2 

J 2 



f y o fit r°° 9 rfr 



< Cy. 

Hence 

G ^ {x > y) - G \V(x,y) + y, y>l. (3 ' 1} 

Let B = (n + 2, oo), n6N. Now assume that n < x < n + 1 and y £ B. We claim that 

G( ,oc)(x,y) < Cx, (3.2) 



where C depends only on a. 
Observe that, 



p6 /*oo 

V(x,y)dy = I p {0<oo) (t,x,y)dydt 



o Jo 



[ P x (r {0jOo) >t)dt<6. 
Jo 



(3.3) 



Consider h(v) = G^ tOQ )(x,v), v G B. By Lemma [2.111 it is regular harmonic on B. Hence 
using the estimate (13. ip we obtain 

G(o,oo)(x,y) = E y G {0:OO )(x,X TB ) 

= E y [G (0 , oo) (x,X TB )-X TB e(0,n + 2 
< E y V(x,X TB ) + C(n + 2). 

Integrating G(p i00 )(v, y) with respect to dv and applying ( 13. 3p we obtain 



n+1 



G i0 , oo) {v,y)dv<6 + C(n + 2). (3.4) 

The final argument for proving (13.21) will use Lemma [2.51 Take D — (n — 1, n + 2), and recall 
that y > n + 2 and x G (n, n + 1) . Due to Lemma 12. Ill the Green function Gto j00 \(u, y) is positive 
regular harmonic on D as a function of u. By Harnack's inequality for harmonic functions on 
D, which follows from Lemma [2.51 we arrive at 

G( ,oc)(x,y) < CG( ,oo){u,y), x,ue (n,n + l), 
which together with (13.4)) completes the proof of the estimate 

G(o,oo){x,v) < Gx, 1 < x <y -2. 
Combining this with (13. ip we obtain 

G(o,oo)(x,y) < C(V(x,y) + x), x > 1. 
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Since G^^(x,y) = x (see (I2.2ip ). then by Lemma 12.101 we have that 

2 

Gr 0)Oa )(x,y) > -x. 

a 

Therefore we proved that 

G(o,oc)(x,y) ~V(x,y) + x, x>l. (3.5) 

To estimate V(x, y) we use 

V(x,y) = f \f>°°\x,y)dt<e & f ' e^°°\x,y)dt 



Jo Jo 

< e 6 G( 0iOo) (x,?/). 

Next consider x < 1 and y < 2. By ( 13.1 ft and Lemma [3.11 we get 

G(p tO0 )(x,y) « G( 0jOo) (z,y) « G( 0>oo) (x, y) + x a/2 . 

Now assume that x < 1 and y > 2. Again G( 0)OO )(-, y), by Lemma 12.111 is regular harmonic on 
(0,2), hence by BHP (see Lemma [23]): 

G(p lOB )(x,y) « G (0;Oo) (l,t/)a; a/2 . 

Due to Theorem 12. 12\ G} Q ^(1, y) < C so by (13.51) we have 

G(o,oo)(l,y) ~ 1, 

which implies 

G(p,oo)(x,y) &x a/2 , x<l,y>2. 
This completes the proof. □ 

Remark 3.3. Lei x < y. Then we have 

{ G {a,oc)( x iy)i x<l, \x-y\< I; 
G \o,oo){ x iV) + x > x>l,\x-y\<l; (3.6) 
x V x a l 2 , \x — y\ > 1. 

4 Exit time properties 

In this section we derive optimal estimates of the expected value of the exit time from a ball 
of arbitrary radius. Then, which seems the most important result of this section, we provide 
optimal estimates of the tail distribution for the exit time from a half-space. That is we improve 
the bounds obtained in Lemmas 12.71 and 12.81 They will play a crucial role in the next section, 
where we deal with the Green function of a half-space in M. d . We start with the one- dimensional 
case. 

Proposition 4.1. For x G (0, R) we have 

E x r im « (x a/2 V x) ({R - x) a/2 V{R- x)) . 
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Proof. If R < 3 then from Lemma 12.51 we have 

E x t^ r) ^(x(R-x)Y' 2 . 

Throughout the rest of the proof we suppose that R > 3. Assume x < R/2. First we prove 
the upper bound. By Theorem 13.21 and (12. 17j) we obtain 

pR i-R 

E x T( ,r) = I G( 0>R )(x,y)dy < J G(p t00 )(x,y)dy 

G\ (x, y)dy + / V /2 V y)dy + (R-x) {x a ' 2 V x) 
Jo 

< 2(aT(a/2))- 1 x a/2 + x(x a/2 Vx) + (R- x)(x a/2 V x) 

< cR(x a/2 Vx). (4.1) 

Now, we deal with the lower bound. By Lemma 12.101 

2 

G(o, R )(x,y) > -G 9 { ^(x,y). 

Denote the first exit time of (0, R) for the Brownian motion by T^ R S y It is well known that 
E x T(q R ^ s = \ x [R — x) (eg. see \±0\ ). Then we have 

-R o rR 



f 2 i 

E x t (0:R) = / G( 10>R )(x,y)dy > - G 9 ^ (x,y)dy 
Jo a Jo 



2 E^- s = -x(R-x). 



a ^ a 
Hence we get, for 1 < x < R/2, 

E x t {0jR) ^xR (4.2) 
Let x < 1. Notice that by the Strong Markov Property 

E x t {0>r) = s(x) + £ a Y (0)2 ), 

where s[x) = E X {E T (°- 2 ) T(o,.r)) is regular harmonic on the interval (0,2) vanishing on its com- 
plement. Therefore by BHP (see Lemma 12.41) we obtain 

s(x) « s(l)x a/2 . 

Moreover due to Lemma [2. 51 we have 

E x r m » x a ' 2 . 

This yields 

E* T{m « (s(l) + l)x a ' 2 . 
Noting that s(l) = E 1 ^^ — E 1 ^^) and observing that (14. 21) implies 

(EV^) - SV (0 ,2)) + 1 « i2, 

we obtain 

^no.fl) ~ i?a;a/2 ' < x < 1. (4.3) 
Putting together (Q|) . fl4T2|) and (@~3]) we obtain 

^(o.fl) « i2 (x Q/2 Vi), for x < i2/2. 
By symmetry we have E x T( 0iR j = E R ~ x T(o )R \, which ends the proof. □ 
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Now we derive bounds for the expected exit times from balls in the multidimensional case. 

Proposition 4.2. For x £ -8(0, R) = {v E M. d : \v \ < R] we have 

E x r mR) « ((R - |x|)"/ 2 V(R-\x\))(RV R a/2 ) . 

Proof. Let be the first exit time from -8(0, R) for the a-stable isotropic process. By the 

result of Getoor [12] we have E x t$q% = c(R 2 — \x\ 2 ) a ^ 2 for c = c(a, d). 
First assume that R < 3. Then from (12.161) we obtain 

E x r B{m « 8^ ( a % = c(i? 2 - |af ) a / 2 « (i?- |x|r/ 2 iT/ 2 , 

which completes the proof in this case. 

Next suppose that R > 3. Let z = x/\x\ if x ^ and 2 = (1, 0, . . . , 0) if x — 0. We now 
take Sr = {v : \ (z,v) \ < R}. The process (z,X t ) is the one-dimensional relativistic process 
(with the same parameter) which starts from \x\. Note that 

E x t B (o,r) < E x t Sr . 

By the one-dimensional result (see Lemma 14. ip we get the upper bound. 

For |x| < R — 1 we get the lower bound by using Lemma \2. 101 and the result for the Brownian 
motion: E x r g B a ^ s R) = ±(R 2 - \x\ 2 ) (see [ID]). Namely 

E x t B {o,r) = I G B{0tR) {x,y)dy>- I G 9 ^L(x,y)dy 
Jb{o,r) a Jb{o,r) 

2 ^(oS) = ^(^-N 2 )- (4-4) 



To complete the proof we need to consider R — 1 < \x\ < R. The conclusion will follow in the 
usual way from BHP (see Lemma [2.4p and the bound above for \x\ = R — 1. We may and do 
assume that x — (0, . . . , 0, \x\). Denote Xq — (0, . . . , 0, R — 1) and zq — (0, . . . , 0, R). Let 

F = -8(0, R) n B(z , 2) 

and 

s(x) =E x E x ^r B(m . 

Observe that s(x) is a positive regular harmonic function on F satisfying the assumptions of 
the second part of Lemma 12.41 hence 

S(x) S3 S(X )(R- \x\) a/2 . 

Next, by the Strong Markov Property 

E x r mR) = s(x) + E x Tp > s(x) + E x t b(x0:1) 

stable 
r B(x ,l) 



s(x )(R~ \x\) a/2 + E x T s E table 



w (s(x ) + l)(R- \x\) a / 2 

= {E x H mR) -E x «T F + l){R-\x\) a l 2 

> cR{R-\x\) a/2 . (4.5) 

The equivalence E x t b ^ u i) ps E x t b ?^ ] s follows from ( 12.16P and 

E xo r mR) - E x °Tp + 1>cR 

follows from f)4.4p . Combining (14.41) and (14. 5p we arrive at the desired lower bound. □ 
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Now, we recall the Ikeda-Watanabe formula [To] which provides a relationship between the 
Green function and the Poisson kernel. Assume D C M. d is a nonempty open set and E is a 
Borel set such that dist(D, E) > 0, then we have 

P' B (X(r D ) (E E,td < oo) = / G?r,(x, - y)^, x E D. (4.6) 

The following generalization of the Ikeda-Watanabe formula was proved in [18J: 

P x (X(T D )eE,t 1 <r D <t 2 )= [ [' p?{x,y)dtv(E-y)dy, (4.7) 

J D Jti 

where < t\ < t 2 , x G D. For D which satisfies the outer cone property we have P X (X TD G 
dD, to < oo) = (see [18J). Therefore the above formulas are true for all sets E C D c for such 
D. In particular, for sets studied in this paper as balls or half-spaces, the process does not hit 
the boundary, when exiting a set. 

As a consequence of formula (14.61) we have the following lemma which proof is omitted. 

Lemma 4.3. Let D <zM. d be a bounded open set then 

Pd{x, z) < E x t d sup u(z-v), z G (D) c , x G D. 
veD 

Moreover, if dist(z, D) > 1 then 

Pd(x,z) < CE x T D e- dist{z ' D) . 
Proposition 4.4. For < x < R we have 

x a / 2 V x 

P X {T{0,R) < T( 0) oo)) « R a/2y R - 

Proof. Assume that R > 1 and < x < R. By Lemma [2. Ill the function G(o,oo)( - , 2i2) is regular 
harmonic on (0, R) , therefore by Remark 13.31 we obtain 

Cx a ' 2 Vx > G (0>oo) (x,2R) = E*G {0)Oo) (X T(0Ry 2R) 

> cE x [X T{OR) A2R;X T(OR) >0] 

> cRP x (X T{OR) >0). (4.8) 

Let n > 3, which we specify later. Again G(o )00 )(-,nR) is regular harmonic on (0,R). 
Applying Remark 13.31 we have 

G( ,oo)(x,nR) = E x G( ,oo)(X noR) ,nR) 

< C (nRP x (X T(0R) > 0)+E*Gl 0>oo) (X T{ORV nR)) . 

(4.9) 
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Moreover Lemma 14.31 and Theorem 12.121 imply 

EXG }o,oo)( X r ((hRV nR) 

G (o,oo) ( v , nR)P {0:R) (x, v)dv 

R 

(n-l)R 

G \o,oo) ( v , nR)P( ,R) (x, v)dv 

R 

oo 



J(n-1)R 

< sup G\ 0tOo) (v : nR)P x (X T(0 G (R, (n - 1)R)) 

R<v<(n-1)R 

poo 

+ sup P {0>R) (x,v) / G} 0oo) (v,nR)dv 

v>(n-l)R J(n-1)R 

< cP x (X T(0R) > 0) + Ce-^ R E x r m) , (4.10) 
where P(o >R )(x, v ) is the Poisson kernel for (0, R) and by Lemma [4.31 it admits 

P(o,R)( x > v ) < CE x r (0;R) e-^ R , v>{n- 1)R. 
Using the (14. 9p and (14. 101) we arrive at 

cnRP x (X T(0R) > 0) > G (0jOo) (x,nR) - Ce~^ R E x r {0tR) . 
By Lemma [4.11 E x T($ tR ) ~ R(x a ^ 2 V x), so Remark 13.31 implies 

G [0>oo) {x,nR) - Ce- {n ~ 2)R E x T {0 , R) > (c - CRe^ n - 2)R ) (x a/2 V 



x 



Now we pick n independently of R > 1 and large enough so that 
c - CRe-^- 2 ^ > c/2. This yields 



„a/2 w 

P x (r {0 , R) < r (0l oo)) > c— — . (4.11) 
Next, for R < 1 we use Lemma [231 to get 

P x (r {0)R) < T(p |00 )) « (4.12) 

Combining (O]) . fliTTD and ff4TT2|) ends the proof. □ 

Now we can prove the main result of this section. 
Theorem 4.5. For x > and t > 1, 

^(r ( o,oc) >*) ~ tl/2 Al. 

Proof. Assume that t > 1. If 2a; > t 1 / 2 the upper bound is trivial. So we may assume that 
2x < t 1/2 . We have 

P x {no,oo) >t)< P x (r (0)R ) >t) + P x (t {0)R) < r (0l oo))- 
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Let R > 2x. By Chebyschev's inequality and Proposition 14.11 we obtain 

P * (TM > t) < £a*s „ 

By the Lemma 14.41 

a;"/ 2 V x 

P X (r^ R) < 7(0,00)) < C Ra /2y R - 

Setting i? = t 1 / 2 we arrive at the upper bound. 
Next, let us observe that by Lemma [2T9l 

/oo /*oo 
p?>°°\x,y)dy> J Eg$$(x,v)dy 

= EP*(T^>T a (t)).° (4.13) 

Let us observe that by Chebyschev's inequality and (12.11) for A = — 1 we have 

P{T a {t) > 2t) = P(e T «W > e 2t ) < Ee T ^e~ 2t = e l ~ 2t < e~ l . 

That is P{T a (t) <2t)>\- e _1 . Taking into account the fact that P^rJ^J > t) « ^ A 1 
we obtain from (I4.13p . 



P *(« >1 »^ E (^P A1 ) 



(()1/2 A1 J Sc (i^ A1 



Now, let x < 1 then 



^(r ( o,oo) > *) > ^r (02) > 0;P X ^^)(r (0iOo) > t) 



> cE 1 



(0,2) 



^(0,2) - Z ' ^1/2 



A 1 



1 



> c(^AljP"(^>2). 



Hence by Lemma 14.41 we obtain 



i^Cno.oo) > *) > 



ex 



a/2 



1 

tV 2 ' 



which completes the proof. 

Corollary 4.6. There exists a constant C such that, for t > and x,y > 1, 

p^\x,y) < C«rW + t-V°) Al)(^Al). 
For t > 1 and x, y > 0, 



□ 



(4.14) 



CrV2 (^ Al ) A l)e-^ <#-)(,,,), 

where c, c\ are some constants. Hence, for x,y,t > 1, satisfying t > \x — y\ 2 we have the 
optimal bound 

#-»(»,») »*-^(^ A l) (^Al' 
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Proof. The upper bound immediately follows from Lemma 12.91 and Theorem 14.51 
Pick < (3 < 1/2 such that (1 + l/(3) a/2 - 2 = 1 and let 

A t = {uo : (3t< T a (t)(u) < 2t}. 

To obtain the lower bound we use again Lemma 12.91 to get 

pf ,C °\x,y)> E g ( °^(x,y);A t 

Next by a classical result 

gl 0,oo \x,y) = 9t{x - y) - 9t{x + y) = 9t{x-y) (i-e~?J 

> ft(x-l/)(lA^)> ft (ar-i/)(^Al)(^Al 

Hence 

pf>°°\x,y) > ct-^e-^P(A t ) (JL Al)(^Al). 
Next we estimate P(Af). By Chebyschev's inequality and by fl2.ll) for A = 1/(3, 



P{T a {t) < (3t) = P( e -W> T ^ > e-*) < e'Ee^^ 



(t) 



e -((l+l//3)«/ 2 -2)t = e -t 



Similarly by (12.1 p for A = — 1, 

P(T a (t) > 2t) = P(e Ta{t) > e 2t ) < e~ 2t Ee Ta{t) = e"*. 

Hence 

P{Al) < 2e~\ 

which implies inf 4> i P{A t ) > 1 — 2e~ 1 and this ends the proof. 



□ 



One of the drawbacks of the inequality in the above Corollary is that the right hand side 
does not depend on the distance \x — y\. The following result will be very useful in the next 
section and it does take into account the distance \x — y\. 

Theorem 4.7. Let x,y > 1 and \x — y\ > 1. Then, for t <\x — y\ 2 , 

P?^ (x, y)<C (r-^-nr A l) (g t (c(x -y)),+tu (c(x - y))) 



\x - y\ 

where c = 8\/2 and C is some constant. Moreover 



r d/2+l,2 V f^\x,y)dt<c{d,a)- 



xy 



\x-y 



d ' 



Proof. Our arguments are based on the idea of proof of Theorem 4.2 in [18]. Throughout the 
whole proof we assume that x, y > 1 and x < y — 1. 
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We first consider the case t < \x — y\ 2 /16. The interval (0, (x + y)/2) we denote by 5* and 
(y — s,y + s) by D(s). Let < s < 1/8, then D(s) G (0, oo)\S. By the Strong Markov Property 
we obtain 

pf' 00 \x,y)dz 

D(s) 

= P%x t eD(s),T {0:Oo) >t) 

< P x (r s <t,X TS >0,X t eD(s)) 

= E x [t s < t,X TS G (0,oo) \ S,P x{Ts) (X t . TS G D(s))] . (4.15) 

Let A = (y - \x - y|/4, y+ \x - y\/4) and B = (0, oo) \ (S U A). Observe that dist(A, S) = 
\x — y\/A and dist(5, D(s)) > \x — y|/8. Because pt{x) is radially decreasing in \x\ we have for 
X TS G B, 



P x ^\X u eD(s)) = [ p u [X TS -z)dz 

Jd{s) 

< \D(s)\ Pu 



lD{s) 

' x - y 



< c|D(s)| ^^j +w (£z| 

where in the last step we applied Lemma 12721 Next observe that gt(x) is an increasing function 
in t on the interval (0, x 2 /2). Hence, for t <\x — y\ 2 /26A, we obtain, for X TS G B, 

P^. ( x t _„ efl(s ))<c| B ( S )|( 9 ,(^) +( ,(|Z|)). 

Define F(t, z) = g t (z / (8^/2)) + tz/(z/(8v / 2))- Then Proposition 14.41 and the above estimate yield 

E x [r s < t,X TS G B,P x ^\X t ^ TS G D(s))} 

< c\D(s)\F(t, x — y)P x (rs < t, X Ts G B) 

< c\D(s)\F(t,x-y)P x (T S <T {0 , oo) ) 

< c\D{s)\F{t,x-y)-^— 

x + y 

< c\D{s)\F(t,x-y)-^-^. (4.16) 

\ x — y\ 

For the set A we have by (14.71) . 

E x [t s < t,X TS G A,P x ^\X t _ TS G D(s))] 

= 11 P?faz) [ u(z -w)P w {X t _ r G D(s))dwdrdz. 

J S JO J A 

Moreover 



P w (X t G D(s))dw = / / p(t,w,z)dzdw 

J a Jd(s) 



D(s) 



pit, w, z)dw j dz 
P w (X t e A)dw <\D(s)\. (4.17) 



D{s) 
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Using (14.171) and observing that v{z — w) < jj((x — y)/4), w G A, z G S we obtain 

E x [ Ts < t,X TS G A,P x ^\X t _ TS G D(s))} 



< c\D{s)\v{{x — y) / A) I I ps(r, x, z)drdz 

Js Jo 



= c\D{s)\v({x-y)jA) / P x (r s >r)dr 

Jo 

< c\D{s)\tv{{x-y)/A) 

< c \D(s)\tj^L_v((x - y)/(8V2)), (4.18) 

where the last step follows from (12. 9p and (12. 6p . Combining (14.151) . (14.161) and (I4.18P after 
dividing by |-D(s)| and passing s \ we obtain for x, y > 1, \x — y\ > 1 and \x — y\ 2 > 256t, 

P?^\x, y) < o^Tyf (<* ((* - 2/)/(8V2)) + ^ ((* - y)/(8V2))) . (4.19) 

Next we consider |x — y\ 2 < 2561 By (14. 14j) we get for t > 1/256 and x, y > 1, 

pf^^^^cr 1 / 2 ^. (4.20) 

Since for t > |x - y| 2 /256 > 1/256, 

cT 1 ' 2 < g t ((x - y)/(8V2)) 

we obtain 

P^°°\x, y) < c°fg t ((x - y)/(8V2)), t > |x - y| 2 /256. 

The above inequality combined with (14.191) , (I4.20p and Lemma 12.21 implies the first claim of the 
theorem. 

To prove the second conclusion of the theorem we apply (14.191) for 256t < \x — y\ 2 and (14.201) 
for 256£ > \x — y\ 2 to get 

> 

r^l 2 pf^{x,y)dt 



< c-^ I [x - yl 7256 r «-w ( Q J x -^y) + tu( x -^S))dt 



\%-y\ 2 A/256 ' \ 9t V Sv^ J 1 V Sv^ 



POO 

+ cxy r d/2 - l dt 

i|z-y| 2 /256 



\x — y\ 2 \x — y\ a / 2+3 \x — y\ c 



xy 

< c- 



\x — y\ d 

Note that we used (12.61) to estimate the density of the Levy measure. □ 
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5 Green function of M C R d , d > 2. 

In this section we extend our one- dimensional estimates for a half-line to higher dimensions. To 
achieve this we start with some upper estimates of the transition densities of the killed process. 
Note that by subordination we have 

p t (x) = Eg Ta(t) (x), xeR d . 

Let A t = {oj : (3t < T a (t)(uj) < 2t} be the set defined in the proof of Corollary 14.61 Let us 
define 

q t (x) = E(g Ta{t) (x);A c t ), x G R d . 
In the sequel we will need a simple upper bound of q(t,x). Note that gt(x) < nj, t > 0, 
and this used for q t (x) yields 

»(*) < ^ P W) < j^e"* (5.1) 

The next lemma will have a very important role in obtaining the upper bound for the Green 
function. We introduce the following notation. For x G M. d we denote x = (x 1; . . . , x^-i) and 
by gt(x) we denote the Brownian semigroup in IR d_1 . 

Lemma 5.1. There is a constant C = C(d,a) such that 

pf(x, y) < Cg 2t (x - y)pi ' oo \x d , y d ) + q t (x -y), x, y G EL (5.2) 

Proof. For y G H and 5 > denote V = V y {5) = [y,y + S] = xf =1 [y h y { + 8} = V x V d C EL 
Then by independence of the subordinator T a (t) and Brownian motion B t one gets 

P x (x t eV,r m >t) 

= P x (X t eV, m > t, A t ) + P x (X t eV, m > t, A c t ) 



E 



x 



A;P X (B To(t) e V|T„(.)) x 



+ P x (X t G V, re > £, 



< sup P X (B U G V)P Xd (B^ (t) G F d , > : < s < t) 



/3t<u<2t 
+ / ^(^ — - 2 ) ( ^- 2 

< CP x (B 2t G V)P^(X, W G ^,r e > t) + / q(t,x- z)dz. 

Jv 

After dividing both sides by \V\ and passing 5 \ we obtain the conclusion. 



□ 



Note that for any x, y G H we can estimate g 2 i(x — y) < ct ^ so from f)5.2p we deduce 
that 

pf(x,y)< C t-^ 2 p ( t °^\x d ,y d ) + q t (x-y), (5.3) 

which will be well estimated with the help of Theorem 14. 71 

Lemma 15.11 the estimate ( 15.11) and Theorem 14.71 show that for the points x, y G EI away 
from the boundary such that \x — y\ > 2 the Green functions for the relativistic process and 
the Brownian motion are comparable. In view of the one-dimensional case this result, proved 
below, is not surprising. 
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Theorem 5.2. For \x — y\ > 2 and x d ,y d > 1 we have 

Proof. The lower bound follows from Lemma [2. 101 
We claim that the following upper bound holds: 

G m (x,y) < c- J,l! '" 



\x-y 



d' 



By gH, 



I^\x d ,y d )<C^., t>l, 



which together with (15.2ft and (15.11) yield the following bound for the transition density 

pf(x, y) < Cg 2t (y - x)|f + ce"> - y\~ d , t > 1. 
Integrating it over (1, oo) we arrive at 

Pt (x, y)dt < C-. -. + | < C Vl 



|y — x|" \x — y\ w \y — x|" 
If \x d — y d \ > 1 we apply (15.31) . (15.11) and Theorem 14. 71 to arrive at 

pf(x,y)dt < C t-«-Wp?>°°\x d ,y d )dt+-^- 

Ji \x-y\ 

< c XdVd 

\%d-yd\ d ' 

Next note that by Lemma 12.21 we can estimate 



pf(x,y)dt< / p t (x - y)dt < 

Jo \x-y 

This combined with ( 15.51) and ( 15. 6ft implies ( 15.4ft . 
Now let d > 3. Since (see (12121) ). 

Gu(x,y) < C- 



d ' 



\x — y 



d-2 



we have the following bound for \x — y\ > 2, 

G M {x,y) < Cminj ^^p , ^J^-a } « Gjf'fay), 

where the last equivalence follows from (12.201) . This completes the proof in this case. 
Now we finish the proof for d = 2. By ( 15.41) . for < 1, we have 

Gu(x,y) < ~ In (l +4^^) = 2nGr ss (x,y), 
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where the last equality is just (12.211) . If > 1> using Lemmas 12.21 and 12.91 together with 

Theorem 14.51 we obtain 



\x-y\ 2 



p t (x - y)dt 

POO 

+ c r l P x {m > t/3)P y (r m > t/3)dt 

J \x—v\ 2 



< 



|z-2/| 2 



< c + C\n 



< Chi 1 + 4 



x — y\ 2 

X2V2 

x — y\ 2 



■ ''■2 !)2 



dt + C r 1 dt + Cx 2 y 2 / t~ 2 dt 



x-y[ 



X2V2 



\x - y\< 



2CirGg tt "(x,y), 



which completes the proof for d = 2. 

Now we are ready to prove the main result of this section. 
Theorem 5.3. For d > 3 and x, y € H; 

\x d V xf)(y d V yf) 1 



Gm{x, y) « min 

G m (x,y) « 
For d = 2 and x, y £ H: 



x d A y d 
k — 2/| 



x — y| d 
0/2 



A 1 



\x — y\ d ~ 2 
1 



\x-y\> 3, 



|x - y\ 



d—a '■ 



x — y\ < 3. 



G m (x,y) w In 1 + 



1 ^x.VxfjfeVj/f) 



|x - 2/| > 3, 



^2 A 2/2 

|x -2/| 



a/2 



A 1 



|x - 2/| 



2-a 



+ ln(l V(x 2 A 2/2)), |x-2/| <3. 



Proof. First assume |x — y| < 3. In the paper [5] it was proved 

G m (x,y) w Gh(x,2/), d>3 

and 

G e (x, 2/) « Gi(x, 2/) + ln(l V (x 2 A 2/2)), d = 2. 
Since |x — 2/| < 3 by Theorem 12. 121 we obtain 

a/2 



Gm(x,y) 



Xd^yd 
\x - 2/1 



and 



Gm(x,y) 



x 2 l\y2 
\x-y\ 



a/2 



A 1 



A 1 



1 



\x - y\ d ~ a ' 



d>3 



\x - y\ 2 - a 



+ ln(l V {x 2 A 2/2)), d = 2. 



□ 
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This yields the bound in the case \x — y\ < 3. 

We introduce the following notation x = (x\, . . . , Xd-i, 1 V Xd), x* = (xi, . . . , Xd-i, 0). Now 
assume \x — y\ > 3 and observe that implies that \x — y\ ~ \x — y\ > 2. 

Then if both points are away from the boundary ( Xd A yd > 1) we use Lemma [5.21 to have 

G m (x,y) « G^ UM {x,y) = G^ uss (x,y). 

Next suppose that x d < 1 < y d . Let = 5(x*, \/2) n BL Then y ^ D(x*) and G H (-,y) is 

a regular harmonic function on D(x*) vanishing on H c . Hence by BHP (see Lemma f2.4j) and 
next by Theorem 15.21 we have 

G m (x,y) « x% 2 Gm(x,y) = x a J 2 G m (x,y) w x^ /2 G^ U88 (x, y). 

A similar argument applies for x^, y<j < 1. Notice that x ^ D(y*) and y ^ Hence 
Gh( - , y) and Ge(^, •) are regular harmonic function on D(x*) and D(y*), respectively, vanishing 
on H c . Hence Lemma [2.41 and Theorem 15.21 imply that 

Gw{x,y) « x a J 2 y^ /2 G m {x,y) « xf 2 y^ /2 G^ uss (x,y). 

Taking into account all cases we have 

G H (x, y) ~ (1 A x d )«/ 2 (l A y d r /2 G 9 m auss (x, y), \x - y\ > 3. 

Applying (12.201) and (12.211) we can rewrite the above bound as 

(xdVl)faVl) 1 ] 
\x — y\ d ' \x — y\ d ~ 2 J ' 

for d > 3, and 

Gh(x, y) « (1 A s^ a (l A y^/ 2 In f 1 + 4 ^1%^ ) , 
for d = 2. Taking into account |x — y| ~ \x — y| > 1, for \x — y\ > 3, we finally arrive at 

I f — y\ \ x ~ y\ I 

and 

□ 

Now we compare the Green functions for half-space for the relativistic process and for the 
corresponding stable process, so we recall the formula of the Green function in the stable case 
(see HI): 

G^ ble (x, y ) = C(a,d)\x-yr d {t + l)d /2 dt - 

One can derive sharp estimates from the above formula. Our results from Theorem l5.3[ Theorem 
13. 21 and Theorem 12.121 show that for the points x, y £ H such that \x— y| < 2 the Green functions 
of the half-space H for the relativistic process and for the corresponding stable process are 
comparable if d > 3. If d — 1 or d — 2 they are also comparable but we have to assume 
additionally that the points are near the boundary. 



G m (x, y )^(lAx d ) a / 2 (lA y d) a/2 



mm 
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Remark 5.4. Suppose that \x — y\ < 2 then 

C G$ able (x,y), d>3; 
G m {x,y)at G^ e (x,y)+ln(lV(x 2 At/ 2 )), d = 2; 

[ G$ able (x,y) + {xAy)V{xA y) a ' 2 , d=l. 

From the estimates obtained in Theorem \5.3\ and Theorem \3. 6 A we can infer that 

G m (x, y) > C(G s ^ able (x, y) + G^ uss (x, y)), x, y G H. 



6 Green functions for intervals 

In this section we provide optimal estimates for Green functions of bounded intervals. We know 
that for any interval the Green function is comparable with the corresponding Green function 
of the symmetric process. That is for the interval (0, R), for R < R , we have 

CiRo^Gftftfay) < G {0 , R) (x,y) < C(R )G s fff(x,y), (6.1) 

where < x,y < R. However, if R grows, then the constant C(Rq) tends to oo, so the above 
bound is not optimal in general case. The aim of this section is to provide optimal bounds for 
large intervals. We recall known estimates for stable cases: 

(5 R (x)S R (y)r/ 2 



111111 { |x-yV Q ' 

ratable f„ „,\ ^ J f -i , (<?fl(^)^fa)) 1/2 



a < 1, 

Gfigfr V)*< ^ (1 + [dR[x Py\ >r " ) > a=1 > ( 6 - 2 ) 

^ mm{(5 R (x)5 R (y)y«-^\ ^fgg^ } , a > 1, 

where 5_r(x) = x A (i? — x). 

We start with the proposition showing that for points x, y in the first half of the interval 
the Green function of the interval and the Green function of (0, oo) are comparable. 

Proposition 6.1. Let R > 4. For x,y < R/2 + 1 we have 

G(o,R)(x,y) ~ G(o iQo) (x,y). 

Proof. Throughout the whole proof we assume that x < y < R/2 + 1. Notice that it is enough 
to prove that 

G {0 ,R)(x,y) > cG(f} t0O )(x,y), 

for x < 1 or | a; — y\ < 1. Indeed, by Lemma 12.101 and Remark 13.31 we obtain for x > 1 and 
\x-y\> 1, 

2 2 1 

G( , R )(x,y) > -G 9 ( ^(x,y) = -x(l-y/R) > —x > cG {0iOo) (x,y). 

We claim that 

G {0iR) (x,y) > CG\ 0)Oo) (x,y), \x-y\<2. (6.3) 
Applying (l2~T4l) with L> 2 = (0, oo) and D 1 = (0, R) we have 



E x 



T(o,R) <T (0 ,oo);e T (°-«'G( 0iOo) (X r(0fl) ,?/) 



< supGjj(z,j/)P x (r (0iK) < T(o )0 o))- 

z>R 
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Next, by Theorem 12.121 and Proposition 14.41 we obtain 

Gfo.cc) (*.!/) " G\ , R) (x,y) < Ce-^R-^ix^ V x). 

Hence Lemma [3.11 yields 

Gf ,co)(^) - G\ , R) (x,y) < Ce~ R / 2 R- 1+a / 2 G\ 0tOQ) (x,y). 

This proves (I6.3P for R > R , if R > 4 is large enough. 

To handle the case 4 < R < R we apply (16.11) together with (16.21) and Theorem 12.121 to 
obtain 

G i0 , R) (x,y) > c{Ra)G%$(x,y) > cG l {0>oo) (x, y), 

which ends the proof of (16. 3p . 

That is, by Theorem 13.21 for \x — y\ < 2, x > 1, we get 

G {0)R) (x,y) > c(Gl oR) (x,y) + G 9 { ^(x,y)) 

> c ( G (o,oo)( x ,y) + x ) > cG i0 ,oo)(x,y). (6.4) 
Next, for x < 1 and y > 2, by BHP (Lemma E3J, 

G( ,R)(x,y) > cG i0)R) {l,y)x a/2 > cx a/2 G { o,oo){^,y) 

> cG(o !OC )(x,y). (6.5) 

Combining H6.4I) and H 6 . 5 1) give us 

G( ,R)(x,y) ~ G(p iOC ) (x, y) , 
for x,y < R/2 + 1. □ 
Proposition 6.2. Lei R > 4. Suppose that l<x<y<R — 1, and \x — y\ > 1 t/jen 

G(o,i?.)(a;, y) ~ G^ } s (x, y) = x(i? - 

Proof. Due to Lemma f2. 101 we only need to prove upper bound. Suppose that 1 < x < y < R—l 
and \x — ?/| > 1. 

At first, let additionally y < 3/4R, then by Remark 13.31 

G( 0jR ){x,y) < G(o,oo){x,y) <cx< Acx(R-y)/R. (6.6) 

By symmetry and the above inequality we have 

G( ,R)(x,y) = G { - Rfi )(-x,-y) = G( ,R)(R- x,R- y) 

< Ac(R-y)x/R, (6.7) 

for R/4 <x< R/2 < y < R - 1 and \x - y\ > 1. 

Hence it remains to consider the case 1 < x < R/4 and 3-R/4 < y < R — l. Denote 
V = r (o,R/2)- Since G^ R )(-,y) is regular harmonic on (0,i?/2), so by Proposition 16.11 and 
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Remark _ 

= E x G(o >R )(X v , y) < E x [X v > R/2; G (0>oo) {R -X V ,R- y) 

< cE x [X v > R/2, \X v -y\<l; G{ 0>oo) (R -X V ,R- y)] 

+ cE x [X„ > R/2; ((R - X v ) V(R- X„) a/2 ) A (i? — y)] 

< c(R - y)P x {r] < r (0>O o)) 

+ cE x [\X V - y\ < 1; G\ 0>oo) (R -X„,R- y)} 

< c{R-y)x/R 

+ cE x [\X V -y\<l; G\ 0>oo) (R -X V ,R- y)] , (6.8) 

where the last inequality is a consequence of Proposition 14.41 Moreover, by Lemma 14.31 we 
obtain 



E x [\X„ - y\ < 1; G\ 0iOo) (R -X V ,R- y)] 
= / P(o,R/2){x,z)G 1 (0oo) {R- z,R-y)dz 

Jy-l 

< cE* V f + e^~ R ^G^ oo) (R-z,R-y)dz 

Jy-l 

p poo 

<c^xe-^- R M J Gl 0iOo) (R-z,R-y)dz 



< cx{R - y)/R, (6.9) 

because y > 3/4R. 

Combining flBTBjl . (IBTjl . (Oj) and fED) we obtain 

G (0) ji)(x, y) < cx(R -y)/R = cG 9 { ^(x, y). 

□ 

Now we can prove the main result of this section. 
Theorem 6.3. Let R > 4 and x < y then we have for \x — y\ < 1, 

G( ,R)(x,y) « mm{G {0:Oo) (x,y),G {0tOo) (R- x,R- y)} 

and for \x — y\ > 1 

(x"/ 2 Vx)((fi-y)°/ 2 V(i?- j /)) 
G(o,R){x,y) « ^ • 

Proof. Observe that by symmetry 

G^ R )(x,y) =G( 0tR )(R-x,R-y). (6.10) 
The y| < 1 follows immediately from Proposition 16.11 Theorem 13.21 and (16.101) . 
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For < x < y < R, \x — y\ > 1 we define x = x V 1 and y = y A (R — 1). Then we can 
repeat the arguments used in the proof of Theorem 15.31 to arrive at 

G {0jR) (x,y) « (lAxr^(lA(R-y))^ 2 Gl^(x,y) 

= (lAx) a / 2 (lA(R-y)) a ^ R ~ y) 

R 

(x a ' 2 Vx)((R-y) a / 2 V (R-y)) , 
= d ' \x-y\>l. 

This completes the proof. □ 
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